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Abstract. We present the results from a detailed deprojection analysis of Abell 1835 as observed by XMM-Newton. 
If we fit the spectra with an isothermal plasma model, the deprojected temperature profile is flat in the outer 
region around 7.6 keV and decreases to ~ 5.6 keV in the center, which may be connected with the gas cooling. 
In the central part, a two-component thermal plasma model can fit the spectrum significantly better. Moreover, 
the cool component (T ~ 1.8 keV) has a much lower metal abundance than the hot component (T ~ 8 keV), 
which may be due to the longer cooling time for the cool gas with lower abundance. In addition, it was found 
that without a main isothermal component, the standard cooling flow model cannot fit the spectrum satisfac- 
torily. From the isothermal model fitting results we also derived the electron density n^, and fitted its radial 
distribution with a double-/? model. The Ue profile inferred with the double-/? model and the deprojected X-ray 
gas temperature profile were then combined to derive the total mass and the total projected mass of the cluster. 
The projected mass is lower than that derived from the weak lensing method. However, assuming that the cluster 
extends to a larger radius ~ 15' as found by Clowe & Schneider (2002), the two results are consistent within the 
error bars. Furthermore, we calculated the projected mass within the radius of ~ 153 kpc implied by the presence 
of a gravitational lensing arc, which is about half of the mass determined from the optical lensing. 
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c/3 1. Introduction the spectrum of the ICM in the cluster center is very help- 

f^ ' ful for exploring the actual physical state and the physical 

^ : There are two yet not very well understood questions con- processes in the center. The second is: whether the cluster 

• --H . cerning clusters of galaxies. The first is: what is the real ^^^^^^ inferred from the two primary observational tech- 

/S ' physical state and process of the X-ray gas in the central ^j^^^g (optical lensing and X-ray observation) are con- 

^ : region of the cluster? The cooling flow rates in the cen- gigtent? For some clusters, the mass determined from the 

■ - - ' ters of most galaxy clusters observed recently by XMM- X-ray method is lower than that from the strong gravita- 

Newton and Chandra are much lower than those predicted ^-^^^^^ j^^gj^g ^jethod by a factor ~ 2 (e.g. A2218, Loeb & 

by the standard cooling flow model (e.g. A1795, Ettori ^^o 1994; A1689, A2163, Miralda-Escude & Babul 1995). 

et al. 2002, Tamura et al. 2001; M87, Matsushita et al. rpj^^ reason for this discrepancy may be a prolate ellip- 

2002). An important implication of the low cooHng flow ^^-^^^^ ^^^^ distribution or a non-isothermal temperature 

rate is that there exist some unknown processes which gtructure existing in the cluster (Miralda-Escude & Babul 

can heat the gas and thus prevent it from cooling. From ^995) ^^ non-thermal pressure supporting the ICM against 

the analysis of the X-ray gas, some authors believe that ^^^^-^^ (^ g^ L^g^^ ^ Mao 1994). Interestingly, when a mul- 

the IntraCluster Medium (ICM) in some clusters is locally ^-^^^^^ analysis is adopted, the X-ray and strong lensing 

isothermal, such as in M87 (except for the regions associ- ^^^gg^g ^^ ^^^ ^o^jij^g g^^ clusters present an excellent 

ated with the radio structures, Molendi 2002), while some agreement, but this discrepancy still exists in some other 

others suggested that the ICM can be better represented clusters which show no evidence for cooling flows (Allen 

by a two-temperature model, in which the cooler com- ^gg^y However, since the standard muhiphase model does 

ponent is associated with the Interstellar Medium (ISM) ^^^ gj^g ^ satisfactory fit to the new Chandra and XMM- 

(Makishima et al. 2001). Apparently, a detailed analysis of Mewton data, it is necessary to explore whether this dis- 

crepancy do exist in the cooling flow clusters. Recently, 
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2 in the cooling flow cluster PKS 0745-191 with XMM- 
Newton observations, while Schmidt et al. (2002) obtained 
a roughly consistent mass for the two methods with the 
Chandra data for Abell 1835. However, on larger spatial 
scales, previous studies have inferred a good agreement 
between the weak lensing mass and the X-ray mass (e.g. 
A2218, A2163, Squires et al 1996, 1997). 

Abell 1835 is an important object to study the above 
properties of the galaxy clusters. It is luminous, with 
a medium redshift (z=0.2523) and a relaxed structure. 
RO S AT a.nd ASCA observations show that it has a large 
cooling flow rate of ~ ^TQOIL^qq MQyr~^ in its cen- 
ter (Allen et al. 1996). However, recent studies based 
on XMM-Newton (Turner et al. 2001) and Chandra 
(Weisskopf et al. 2000) show that the mass deposition rate 
is not that large. The RGS onboard XMM-Newton has lim- 
ited the cooling flow rate to 315 M0yr~^ within a radius of 
150 kpc (Peterson et al. 2001), and Chandra derived that 
the cooling flow rate is ^ 500 M0yr~^ within a radius of 
250 kpc (Schmidt et al. 2001). Furthermore, the XMM- 
Newton RGS did not detect any multiphase gas below a 
certain low temperature {'^ 2.7 keV), inconsistent with the 
prediction of the standard cooling flow model (Peterson et 
al. 2001). The gas temperature of Abell 1835 was deter- 
mined as ~ 12 keV from the Chandra data (Schmidt et al. 
2001), but was determined as ~ 7.6 keV from the XMM- 
Newton data (Majerowicz et al. 2002). The discrepancy 
between the model prediction and observations as well as 
that between different observations urges a more careful 
examination with high quality data. 

We investigate here the temperature, density, cooling 
flow rate and mass of Abell 1835 with a deprojection tech- 
nique based on the data observed by XMM-Newton EPIC. 
The deprojection technique can reveal the real spectra of 
the cluster gas in different spherical shells, and we can 
further determine the deprojected temperature and the 
mass distribution of the gas in the cluster (e.g. Chen et 
al. 2003). The XMM-Newton EPIC is the most sensitive 
X-ray telescope which also has high spatial and spectral 
resolutions, and therefore meets all the requirements for 
the detailed spectral analysis. 

The structure of this paper is as follows: Sect. 2 de- 
scribes the observation, background correction and the 
spectral deprojection technique. Sect. 3 presents the de- 
projected spectral analysis with three different models: 
single-temperature model, two-temperature model and 
cooling flow model. In Sect. 4 we obtain the electron den- 
sity profile, calculate the total mass, total projected mass 
and the projected mass within the optical lensing arc, then 
we discuss the discrepancy between the X-ray mass and 
the optical lensing mass. Our conclusion is given in Sect. 
5. 

Throughout this paper, the energy band is 0.5 ^ 10 
keV, and unless otherwise noted we use a cosmology with 
iJo = 50 km s"iMpc"\ go = 0.5, VL^ = 1.0, and f^A = 0. 
Therefore, 1' corresponds to 296.6 kpc at the distance of 
AbeU 1835. 



2. Observation and data preparation 

Abell 1835 was observed by XMM-Newton during the 
phase of performance verification (observation ID is 
0098010101). Since MOSl was operating in Large Window 
mode, we only use the data coming from M0S2 and pn 
cameras which were operating in the Full Frame mode and 
with the thinl filter. The total exposure time is 60 ksec, 
of which only about 26 ksec is usable for spectroscopic 
analysis for both pn and M0S2. For the M0S2 data, we 
use the event with PATTERN < 12, and for the pn data 
PATTERN < 4. The calibration is performed in SAS 5.3.3. 

2.1. Background correction 

Background subtraction was carried out using the same 
method as that of Majerowicz et al. (2002). The pn data 
have another source of contamination called out-of-time 
(OCT) events counted during the read-out (see Striider 
et al. 2001). This contamination has also been corrected 
for in our analysis. 

2.2. Spectral deprojection 

Abell 1835 appears to be a relaxed cluster of galaxies, 
therefore we assume that the temperature structure of 
this cluster is spherically symmetric, and apply a depro- 
jection technique as done by Nulsen & Bohringer (1995). 
We divide the image of the cluster into 7 annular regions 
centered on the emission peak for the extraction of spec- 
tra and use the outmost ring (8.33' — 10.42') to determine 
the local Cosmic X-ray Background (CXB). However, in 
the sixth ring the signal-to-noise ratio is low, so we only 
consider the inner five regions (r < 6'). The minimum 
width of the rings was set to 0.75' which is wide enough 
to ignore the PSF (Point Spread Function), whose FWHM 
(FuU Width at Half Maximum) is 5" for M0S2 and 6" for 
pn. For each annular region, an Ancillary Response File 
(ARE) is generated using SAS, then through the ARE, the 
vignetting (Arnaud et al. 2001) correction is administered. 
Deprojected spectra are calculated by subtracting the 
contribution from the outer regions for all spectral com- 
ponents. Within each annular region, the spectrum per 
unit volume is assumed to be the same. The deprojected 
spectrum of the ith shell is then calculated by subtracting 
the contributions from i+lth to the outmost shell from 
the annular spectrum of the corresponding radius (e.g. 
Matsushita et al. 2002). 



3. Spectral analysis 

3.1. Single temperature model 

We analysed the deprojected spectra of both M0S2 and 
pn data using XSPEC version 11.2.0 (Arnaud 1996), and 
we selected the following model: 



Modeli = Wabslnn) x Mekal{T,z,A,norm), 



(1) 
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here Wabs is a photoelectric absorption model (Morrisson 
& McCamraon 1983) and Mekal is a single temperature 
plasma emission model (Mewe et al.l985, 1986; Kaastra 
1992; Liedahl et al. 1995). We fixed the redshift z to 
0.2523 and nn to the Galactic absorption 2.24x10^" cm^^ 
(Dickey & Lockman 1990). The results are listed in Table 
1 and the spectrum of the central region (r < 0.75') for 
M0S2 data fitted by this single temperature model is 
shown as Fig. 3 (a). 

To derive the temperature profile of this cluster, we 
fitted the combined spectra of M0S2 and pn with a single 
temperature model; the results are also listed in Table 1. 
From the deprojected temperature profile (the diamonds 
in Fig.l), we can see that the temperature is nearly con- 
stant in the outer regions but decreases towards the center, 
which may be connected with the gas cooling. Then we fit 
the temperature profile with the formula: 



T{r) = 



Jjr 



(2) 



The best fit parameters are: a — -101.9 keV, h — -0.11 
arcsec"^, c = 7.55 keV. 

Our result does not differ much from that of 
Majerowicz et al. (2002) based on the same data, al- 
though what Majerowicz et al. derived is the projected 
temperature. However, the deprojected temperature from 
the Chandra data (Schmidt et al. 2001), ~ 12 keV, is much 
higher than our result. This may be due to the difficulty 
in identifying the background flares in the Chandra data 
(Markevitch 2002). 
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Fig. 1. Temperature profile of Abell 1835 for the combination 
of MOS2 and pn with a confidence level of 90%. (diamonds: tem- 
perature from the IT model fitting of the combined spectra; 
stars: temperature from the 2T model fitting of the pn central 
spectrum; triangles: temperature from the 2T model fitting of 
the M0S2 central spectrum. We have offset the triangles 10" 
to the left so as to illustrate the two fits clearly.) The solid line 
is the best-fitted profile of the IT model. 



3.2. Double temperature model 

From the null hypothesis probability P (the adopted level 
of a fit is P > 0.1, and when 0.01 < P < 0.1, the fit is 
marginally acceptable) in Table 1 we find that the spectral 
fitting with the single-temperature model is acceptable for 
the outer regions, but unacceptable for the central region. 
Since the calibrations of pn and M0S2 are not quite con- 
sistent with each other, we fit the central spectrum of pn 
and M0S2 with the two-temperature model respectively: 

Model2 — Wabs^nn) X {Mekal{Ti, z, Ai,normi) + 

Mekal{T2, z, A2,norm2)), (3) 

which means that there are two components in this region 
with different temperatures, and the free parameters have 
the same meaning as in Sect. 3.1. Assuming that the two 
components have the same abundances, we fix the abun- 
dances of the cool component to those of the hot gas (see 
Fig. 3 (b) for M0S2), and obtain the best-fit parameters 
listed in Table 2. Note that the cool component only has 
a small volume filling fraction of less than a few percent 
for both pn and M0S2 (see f^oi in Table 2). Although the 
best parameters derived from pn and M0S2 do not agree 
with each other in detail, the values of x^ for both of them 
improved significantly and are marginally acceptable. 

We further leave the abundances free and obtain a even 
better fit (see Fig. 3 (c)). The resultant model parameters 
of pn and M0S2 arc now consistent with each other, as 
shown in Table 3, where F is derived from a simple _F-test 
(Bcvington 1969) and indicates the degree of confidence of 
the model in which both abundances are left free param- 
eters compared to the model setting the two abundances 
the same. From F we know that the substitution of this 
model for pn data is significant with > 99.9% confidence, 
while for M0S2 it is about 80%. The existence of the two 
abundance components implies that the distribution of 
the abundance for both temperature components is inho- 
mogeneous in the central region of the cluster. As shown 
in Table 3, in the steady cooling fiow scenario, there ex- 
ist two temperature components, and the cool one with 
a temperature less than 2 keV has a lower abundance. 
This may be due to the metal-rich gas cooling together 
with the metal-poor gas by bremsstrahlung emission until 
^2 keV, below which line cooling becomes important, and 
the metal-rich gas would cool at a much faster rate than 
the metal-poor gas (Fabian et al. 2001). So in the steady 
cooling flow scenario the metal-rich component cannot be 
seen at temperatures below about 2 keV and thus the cool 
component presents a lower abundance. 

Although the combined spectra of pn and M0S2 are 
not suitable for the two-temperature model fitting because 
of the slight inconsistence between their calibrations, we 
still perforin the combined fits to calculate the mass in 
the 2T model method easily (see Sect. 4.4). The result- 
ing parameters as listed in Table 2 and 3 lie between the 
two sets of parameters derived separately from the pn and 
M0S2 data. 
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Table 1. The best-fit free parameters of Abell 1835: the temperature T of M0S2, pn and the combination of them; the 
abundance A and the normalized constant norm for the combination of MOS2 and pn. norm — 10^^"'' / {4-kD^) j nenndV , 
where D is the distance to the source (cm) and n^ is the electron density (cm""^). Lx is the bolometric luminosity and P is the 
null hypothesis probability of the spectra fitted in XSPEC. The errors represent a confidence level of 90%. 
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Table 2. The best-fit parameters of Abell 1835 using a two-temperature model for the central region(0.0' — 0.75') and assuming 
that the two components have the same abundance. The error bars are at the 90% confidence level, fvoi is the volume fraction 
of the cool component. 



Detector 



T (kcV) 
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(solar) normi norm2 
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Table 3. The best-fit parameters of Abell 1835 using a two-temperature model for the central region (0.0' — 0.75'), but leaving 
the abundances of the two components free. The error bars are at the 90% confidence level. F is the confidence level of using 
this 2T2A model instead of the 2T1A model. 
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We also fit the outer spectra with the two-temperature 
model, but the fits do not improve significantly. 

From the analyses above, we know that a two- 
temperature model may represent the spectrum of the 
central region much better, which is possibly due to 
the existence of cooling gas in the cluster center or the 
possible presence of ISM associated with the cD galaxy 
(Makishima et al. 2001). However, for the outer parts, 
the fit of the single-temperature model is good enough 
to be acceptable. So we inferred that in the cluster cen- 
ter there exist two gas components with different temper- 
atures, while in the outer regions the gas is isothermal 
(Kaastraet al. 2003). 

Fig. 2 is the deprojected abundance profile. It can be 
seen that the abundance is higher in the cluster center 
and tends to constant in the outer region. This indicates 
that the excess metal in the cluster center is produced in 



the cD galaxy and ejected into the ICM (Makishima et al. 
2001). 



3.3. Asorbing column density 

To verify if the observed nn is consistent with the Galactic 
column density, we fit the spectra with the parameter nn 
left free. First, we perform joint fits to 5 pn and M0S2 
spectra respectively with a single temperature model, as- 
suming that the 5 spectra have the same uh, and the 
results are shown in the upper part of Table 4. It can be 
seen that the two best-fit nj^ are much smaller than the 
Galactic column density ~ 2.24 x 10^° cm^^. Since the 
central spectrum can be fitted well by a two-temperature 
model, we consider a joint fit with a two-temperature 
model for the central spectrum and a single temperature 
model for the outer 4 spectra, and the hydrogen column 



Jia et al.: The Analysis of Abell 1835 



1.U 








- 


0.8 


- 






- 


0.6 


- 








- 




; 


S( 






- 


0.4 


- 


A 






- 




Y 
















.. 


_ 


0.2 


- 






' 








/ 


^ 
i 


i 








- 


00 













50 100 

Radius (orcsec) 



150 



Fig. 2. Abundance profile of Abell 1835 in units of solar metal- 
licity with a confidence level of 90%. The symbols have the 
same meaning as Fig.l. 



density nn of each spectrum is assumed to be the same. 
The results are listed in the lower part of Table 4, where 
F is the confidence level corresponding to using the 2T 
model to fit the central spectrum instead of using the IT 
model to fit all five spectra. We find that the fits become 
better, and from F we know that the substitution by a 2T 
model is necessary. It is also shown in Table 4 that the hh 
determined from the pn spectra is still smaller and about 
half the Galactic column density, while the nn obtained 
from the M0S2 spectra is consistent with the Galactic col- 
umn density. We should note that for many clusters the 
best-fit column density derived from the MOS is always 
in good agreement with the Galactic value, while that es- 
timated from pn data is systematically smaller (Molendi 
& Pizzolato 2001). This may be due to the MOS detector 
having a more reliable calibration than the pn detector. 

We also note here that the joint fit with a two- 
temperature model for the central spectrum can yield a 
much more believable nn, and this is also evidence for 
the existence of two gas components with different tem- 
peratures in the cluster center. 

Through the analysis above, we can conclude that 
the observed uh of XMM-Newton is consistent with the 
Galactic column density, and so the analysis in Sect. 3.1 
and Sect. 3.2 is reliable. In addition, the MOS data are 
probably much better for determining nn than the pn 
data. 



3.4. Cooling flow model 

X-ray observations of clusters of galaxies show that in the 
central regions of some clusters the cooling time of the 
ICM is significantly less than the Hubble time (Edge et al 
1992; White, Jones & Forman 1997; Peres et al. 1998). 
ROSAT and ASCA observations show that AbeU 1835 
is the cluster that contains the largest cooling flow rate 



Table 4. The absorbing column density determined by pn and 
MOS2 spectra assuming that all the spectra have the same 
riH- IT model: the joint fit with a single temperature model. 
2T-I-1T model: the joint fit with a two-temperature model for 
the central spectrum and a single temperature model for the 
outer 4 spectra, where F is the confidence level of the substi- 
tution by the 2T-I-1T model. 



Model Detector 
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(10^ 



Xved/dof 



F 
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pn 
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2T+1T pn 
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01+"-^" 
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1.97Z 



1.15/1114 0.0004 
1.15/385 0.02 

1.06/1111 0.08 >0.999 
1.13/382 0.035 >0.90 



(Allen et al. 1996). Here we use the M0S2 data, which 
are assumed to be better for the determination of nn and 
thus the cooling flow rate (Molendi & Pizzolato 2001), 
to calculate the cooling flow rate of Abell 1835. There are 
two different methods for calculating the cooling flow rate: 
spectral method and spatial method. 

Since we have already deprojected the spectrum, i.e., 
the contribution of the ambient gas has been excluded, the 
spectrum should be fitted well with the standard cooling 
flow model: 

Modeh = Wabsinn) x {ZwabsiAnn) x Mkcflow{M))iA) 

where Wabs was described in Sect. 3.1, Zwabs is an 
intrinsic photoelectric absorption model (Morrison & 
McCammon 1983), and Mkcflow is a cooling flow model 
(Fabian, 1988); Auh is the intrinsic absorption and M the 
rate of gas cooling out of the flow. Since the cooling flow 
radius of Abell 1835 is - 230 kpc (Allen 2000) at which 
the cooling time first exceeds the age of the universe, we 
only fit the central spectrum (~ 225 kpc) with this model 
(Table 5) . It can be seen that to fit the spectrum a cut-off 
temperature, lowT, is needed, which is found to be 2.3 
keV. If we fix lowT — 0.01 keV, the x^ value will be much 
larger (see Table 5) and from P we know that this fit is 
unacceptable. Therefore, the standard cooling flow model 
cannot fit the spectrum satisfactorily. 

When adding an isothermal Mekal component, the fit 
becomes much better, as shown in Table 6. F shows that 
the replacement ((Mekal-|-Mkcfiow) instead of Mkcflow) 
is significant, i.e., the cooling flow model with a Mekal 
component is more reliable, which implies that there 
should exist an isothermal component in the cluster. Then, 
the coofing flow rate M is about 656.2t3go;2 Moyr^i 
within 225 kpc (0.75'). The spectrum of the central re- 
gion (0.0' - 0.75') fitted by this model for M0S2 data is 
shown in Fig. 3 (d). 

It can also be seen from Table 6 that the intrinsic ab- 
sorption Ann of Abell 1835 is very small and consistent 
with zero within the error bars, which is in agreement 
with the previous results, such as those of RGS (Peterson 
et al. 2001) and Molendi & Pizzolato (2001) in which the 
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Table 5. The best-fit parameters for the central region (r < 0.75') of M0S2 data by the standard cooling flow model. The 
errors are at the 90% confidence level. M is the mass deposition rate, and Ann is the intrinsic absorption. 
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Table 6. The same as Table 5 but for the cooling flow model with a Mekal component. F is the confldence level when using 
(Mekal±Mkcflow) instead of the standard cooling flow model. 
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(keV) 
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(10 



norm 

3„„- 



-o.aii 

-1.13 
-0.63 
-0.85 



M 

(Mo; 



+403.4 
360.2 
169.5 



Ahh 
(10^^ cm-^) 



Xvcd/ do f 



F 



7.30 
6.38 



0.01 (fix) =r„ 

0.01 (fix) =Trr. 



0.44 ± 0.08 
0.43 ± 0.07 



7.20; 
8.4i: 



656.2 
262.5 



0.06 ± 0.06 
0.25 (fix) 



1.22/152 
1.31/153 



0.03 
0.006 



>0.999 



projected spectra are used. If we adopt the best-fit result 
from the analysis in Schmidt et al. (2001) of the Chandra 
dataset, and fix Ann = 0.25 x 10^^ cm~^, we measure 
M = 262.5^262:5 M0yi'"\ which is unacceptable at the 
adopted level of P=0.01. So the intrinsic absorption of 
Abell 1835 is very small and close to zero. 

Another method for calculating the cooling flow rate is 
the spatial method. According to the energy conservation, 
this method can be expressed as (White, Jones & Forman 
1997): 



J2 M{i')[Ah{i) + Acj,{i)], 



(5) 



where L^ is the bolometric luminosity, shown in Table 
1; M{i) is the mass deposition in shell i; X]i'=i~ M{i') 
is the mass of gas that needs to pass through shell i to 
give rise to the radiation and mass deposition in interior 
shells; A(f){i) is the change in the gravitational potential; 
h{i) — ^kT{i)/ iinip is the temperature in units of energy 
per particle mass of the hot gas; f{i) is a fraction of the 
overall change in the cluster potential, and in this analysis 
f(i) = 1 is used. The first two terms in this equation 
represent the mass that is left in shell «, and the second 
two terms represent the mass that flows through to the 
next interior shell. From this method, we can estimate the 
cooling flow rate within the central region (r < 225 kpc 
^ about 1600 Moyr^^, consistent with the previous result 
of Chandra from the same method (Schmidt et al. 2001). 
The result derived from the spatial method gives the 
total loss rate of the energy, but the cooling flow rate 
obtained from the spectral method only contains the mass 
deposition rate during the cooling flow. The discrepancy 
between these two methods implies that there should exist 
some other energy sources, for example heating by AGN 
(see e.g. Bohringer et al. 2002) or thermal conduction (e.g. 
Voigt et al. 2002) to support the X-ray emission in the 
center of the cluster. 



With any acceptable model ( the two-temperature 
model or the cooling flow model with a mekal compo- 
nent) we have investigated above, we can conclude that 
in the cluster center there should always exist a locally 
isothermal component and there should also exist a minor 
component which may be a single phased cool component 
or a multiphased one. However, with the current data we 
can not discriminate between these two models. 

4. Mass analysis 

4.1. Electron density 

Here, we have divided the cluster into 17 annular regions 
centered on the emission peak. By the deprojection tech- 
nique we can calculate the photon counts in each corre- 
sponding shell. Then, using the abundance and the tem- 
perature profile, the normalized constant of each region, 
norm, can be estimated. Furthermore, since: 



norm = lO^'^'^/iAnD^) / n^nndV, 



(6) 



we can derive the deprojected electron density ng, shown 
as the triangles in Fig. 4. 

We fit the electron density with both single-/? and 
double-/3 model. First, we fit it with a single-/? model (e.g. 
Cavaliere & Fusco-Femiano 1976) in which the electron 
density profile ne{r) is defined as: 



) ne{r) = no 



1± - 



(7) 



where rig is the central density, Tc is the core radius, and (3 
is the slope. The best-fit parameters are: no = 0.06 cm^^, 
re = 0.21 arcmin, /3 = 0.56, x^=60.33, D0F=14 and the 
best-fit single-/3 model profile is plotted as a dotted line in 
Fig. 4. We can see that it docs not fit well, especially in the 
outer regions. For this reason we fit it with a 'double-/?' 
model (see Chen et al. 2003): 



ne{r) = rioi 



1+1 — 
rci 



ft 



-no2 



r 
rc2 



ft 



.(8) 



Deprojected Spectrum; NH=0.224E21; 1T 
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Deprojected Spectrum; NH=0.224E21; 2T1A 




channel energy (keV) 
Deprojected Spectrum; NH=0.224E21; 2T2A. 



ctionnel energy (keV) 
Deprojected Spectrum; NH=0.224E21; CoolingFlow 
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Fig. 3. The spectrum of the central region (r < 0.75') for M0S2 data of Abell 1835. a) fitted by the IT model; b) fitted by the 
2T1A model; c) fitted by the 2T2A model; d) fitted by the cooling fiow model (wabs(mekal+zwabs(mkcflow))). 



Note that this model fits the electron density directly and 
is different from the common double-/? model that was 
used to fit the surface brightness. The best fit parameters 
are: noi = 0.003 cm~^, Vd — 1.84 arcmin, /3i = 1.24, no2 
= 0.08 cm-3, rc2 = 0.155 arcmin, P2 = 0.60, x^ = 12.25, 
D0F=11 and the best- fit line is shown as the solid line in 
Fig. 4. It can be seen from the x^ value that the electron 
density profile is fitted better by a double-/? model than 
by a single-/? model, so in the following calculations we 
use the electron density profile inferred from the double-/? 
model. 

4.2. Total mass 

Assuming spherical symmetry and hydrostatic equilib- 
rium, the total mass profile can be determined when the 
radial profiles of the gas density and temperature are 
known. We calculate the gravitational mass of Abell 1835 
with the hydrostatic equation (Fabricant et al. 1980): 



Mtoti< r) = - 



fcsTr^ d(lnne) rf(lnT) 
Gimip dr dr 



(9) 



here ks is the Boltzmann constant, G is the gravitational 
constant, /x is the mean molecular weight of the gas in 




0,1 1,0 10-0 

Radius (arcmin) 

Fig. 4. Electron density profile of Abell 1835. The error bars 
are at a 68% confidence level. The solid line is the best-fit 
profile for a double-/? model, and the dotted line for a single-/? 
model. 



units of nip, the proton mass. For a fully ionized gas with 
a standard cosmic abundance, a suitable value is /^ = 0.6. 



Jia et al.: The Analysis of Abell 1835 



Using the best-fit profiles of tlie electron density ne{r) 
and the deprojected temperature T{r) we can obtain the 
total mass profile shown in Fig. 5. It was found that the 
total mass Mtot within the radius of 6' is 1.05±0.13 x 10^^ 
Mq, which is consistent with that found by Majerowicz et 
al. (2002) within the errors. 



E 

3 0-4 
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Fig. 5. The total mass profile of Abell 1835. The error bars 
(dotted lines) are at a 68% confidence level. 
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Fig. 6. The projected mass profile in the universe model: Ho 
= 100 kms"^Mpc"\ fim = 0.3, Qa = 0.7 (shown as the sohd 
line; the dotted lines mark the confidence level of 68%), and 
compared with that of Clowe & Schneider (shown as the dashed 
line). In the upper panel, the solid line presents the projected 
mass within 6'; the lower panel is for the projected mass within 
15'. 



4.3. Total projected mass 

We calculate the projected mass profile of Abell 1835 for 
another universe model: Hq — 100 km s~^Mpc^^, 51„i ~ 
0.3, and JIa = 0.7 and compare it with that of Clowe & 
Schneider (2002), who also calculated the total mass in the 
same universe model using the weak lensing method based 
on the data from the ESO/MPG Wide Field Imager (see 
Fig. 6). The upper panel shows our projected mass within 
6', and it can be seen that the two results are consistent 
within our error bars in the radial range between 1' and 
4', where lensing mass estimates are available. By increas- 
ing the radius, our mass measurement becomes lower than 
that of Clowe & Schneider. Since Clowe & Schneider have 
derived the lensing mass profile up to 15', we also extrapo- 
late our best-fit temperature and electron density profiles 
to that large radius, then calculate the projected mass, 
shown in the lower panel. The two results are consistent 
within our error bars. Moreover, if we assume that the 
mass in the outer regions is distributed as derived from 
Clowe & Schneider, we estimate a contribution projected 
from r > 6' to the inner part of about 0.16 x 10^^ Mq. By 
considering this contribution, the two results in the up- 
per panel will also be consistent. Therefore we conclude 
that the X-ray mass shows an excellent agreement with 
the weak lensing mass at large radii. 



4.4. Mass within the optical lensing arc 

The optical observation (taken on 27 February 1998 with 
the Canada France Hawaii Telescope, CFHT) shows an 
arc ~ 153 kpc from the center of Abell 1835 which is 
thought to be a gravitationally lensed image of a galaxy far 
behind AbeU 1835 (see Fig.2 of Schmidt et al. 2001). We 
estimate the X-ray projected mass inside this arc: Marc = 
0.80 ±0.02 X 10" Mq (errors are at 68% confidence level). 

Since there is a two-temperature structure in the cen- 
tral part, we also calculate the projected mass in the 2T 
model: assuming that the two temperature components 
are in pressure equilibrium. Through the two-temperature 
model fitting, we obtained the temperature and norm for 
each component. Then using Eq. (6) the electron density 
Tie of each can be derived. Because of the pressure equilib- 
rium, we can calculate the projected mass from any com- 
ponent by using Eq. (9), and the result is about 0.7 x 10^^ 
Mq, which is still in agreement with that from the IT 
model. 

Allen et al. (1996) calculated the lensing mass in- 
side the radius of this arc, and the result is between 
1.4x10" Mq and 2x10" Mq for the redshift of the arc 
0.6 < Zarc < 3.0. Schmidt et al. (2001) also calculate the 
X-ray projected mass inside this arc using the Chandra 
data and their value is 1.06^1.27x10^^ Mq. Both results 
are larger than ours. We will discuss the reasons in the 
following section. 



Jia et al.: The Analysis of Abell 1835 



4.5. The discrepancy between the X-ray mass and the 
optical lensing mass 

Sect. 4.4 shows that the projected mass within the opti- 
cal lensing arc we calculated is not consistent with previ- 
ous results. Our estimate is lower than that derived from 
the analysis of Chandra data in Schmidt et al. (2001), 
the main reason being probably the large difference of the 
temperatures determined by XMM-Newton and Chandra, 
which directly causes the difference of the mass. The X- 
ray gas temperature determined by Chandra is about 12 
keV, but it is only about 7.6 keV for XMM-Newton. 

On the other hand, Bartelmann (1995) shows that on 
average the cluster mass required for large arcs will be 
lower by a factor of ^ 1.6 than expected from radially 
symmetric models, and that there is a probability of ^ 
20% of overestimating the actual cluster mass by a factor 
of ~ 2. The inner part of the X-ray isocontours of Abell 
1835 is elliptical with an axis ratio of / = 0.85 (see Fig. 4 
in Schmidt et al. 2001). Therefore, one reason for the dif- 
ference between our result and the lensing mass might be 
the assumption of a symmetric model when calculating 
the lensing mass. Another reason may be the existence of 
radio structure (Ivison et al. 2000) in the center region of 
Abell 1835 like in PKS 0745-191 (Chen et al. 2003), al- 
though the radio emission of the latter is much stronger. 
The radio plasma may fill a larger volume in the central 
region of the cluster or provide some additional pressure 
to support the X-ray gas. Therefore, when we consider the 
radio component, the projected mass may be larger and 
may be consistent with that of the optical lensing. 

5. Conclusion 

We have presented a detailed analysis of about 26 ksec 
of XMM-Newton observations of the galaxy cluster Abell 
1835. Through the deprojected spectra analysis we derived 
the deprojected temperature profile, which is flat in the 
outer regions but decreases towards the center. The abun- 
dance is higher in the cluster center, which may be caused 
by the cD galaxy (Makishima et al. 2001). We also fit the 
spectrum in the central region with a two-temperature 
model, and the fit becomes better. It was found that the 
two temperature components coexist in the cluster cen- 
ter which is due to the gas cooling and/or the ISM asso- 
ciated with the cD galaxy (Makishima et al. 2001). The 
lower temperature component also has a lower abundance, 
which may be due to the fact that the metal-rich gas cools 
faster when temperature falls below about 2 keV (Fabian 
et al. 2001). Then considering a cooling flow model, we find 
that the intrinsic absorption of Abell 1835 is very small if 
not zero and the cooling fiow rate is 656.2l3gQ'2 MQyr~^, 
which is much smaller than that derived by the spatial 
method of ^ 1600 Mq yr~^. This implies that there should 
exist some other energy sources that can heat the gas and 
prevent it from cooling down. Furthermore, we showed 
that the cooling flow model without a Mekal component 
cannot explain the data satisfactorily. Therefore, with any 



acceptable model which we investigated in Sect. 3, the 
main component in the cluster center (r < 0.75') is always 
a locally isothermal component and there should also exist 
another minor component, which may be a single phased 
cool component or a multiphased one. However, with the 
current data we cannot discriminate between these two 
models. 

We fit the electron density profile by a double-/? model, 
and use the best fitting parameters to calculate the to- 
tal mass. We find that within a radius of 6', Mtot = 
1.05±0.13x 10^'"' Mq, which is consistent with the value 
derived by Majerowicz et al. (2002) within the error bars. 
We also calculate the total projected mass using a Hq 
— 100 km s^^Mpc^^, iljn — 0.3, fl\ — 0.7 cosmological 
model and compare it with what was calculated in Clowe 
& Schneider (2002) by the weak lensing method. We find 
that the X-ray projected mass within 6' is slightly lower 
than the weak lensing mass. However, if we assume that 
our best-fit ne{r) and T{r) can be extrapolated out to 
15' where mass measurements from Clowe & Schneider 
(2002) are available, we show that the two determinations 
are consistent within the error bars. 

Lastly, we calculate the projected mass within the op- 
tical lensing arc (r ^ 153 kpc) using the single tempera- 
ture model, which is 0.80 ± 0.02 x 10^^ Mq and is almost 
the same as that from the double temperature method. 
This value is about half of the gravitational lensing mass. 
Our result is similar to that of another large cooling flow 
cluster PKS 0745-191 (Chen et al. 2003). Therefore, our 
results imply that the mass discrepancy between the X- 
ray and the optical lensing does exist in the well relaxed 
clusters. This discrepancy may be due to (i) the assumed 
symmetry adopted in calculating the lensing mass, and 
(ii) a pressure support by the radio plasma in the cluster 
center. 
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